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A NEW ALGORITHM FOR ATTITUDE-INDEPENDENT 
MAGNETOMETER CALIBRATION 


Roberto Alonso* and Malcolm D. Shuster t 


A new algorithm is developed for inflight magnetometer bias deter- 
mination without knowledge of the attitude. This algorithm com- 
bines the fast convergence of an heuristic algorithm currently in use 
wuh the correct treatment of the statistics and without discarding 
data. The algorithm performance is examined using simulated data 
and compared with previous algorithms. 


Introduction 


At orbit injection, the only attitude sensor which may be operating is often the vector maa 
i* c “ in r* s P acecraft - Thus, some means is usually needed to quickly determine thk 

STdf tad +2SZZTT ° f ,he * “”" 0 ' atSage 

me des,red a| g° nthm must not require a knowledge of the attitude as input. 

A number of algorithms have been proposed for estimating the magnetometer bias The 
“7'“! ? 10 “** f?I N* vector by minimizing ute weight suZCTuares otreS 

magKhfflelds M t ' ThLT “ ° f ma e" itudes oflhc measure and modeled 

.. ® ds 1 1 ] • This approach has the disadvantage that the cost function is nnartie in 

the magnetometer bias, and therefore admits multiple minima. If these soludZ are cl^ o 

ne another, t hen convergence of the algorithm may be poor. Typically, one initiates the least- 
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So— 

field. 

r^mhhir r 1 21 advocated centering the data to remove the quartic dependence. This leads 
- TSM t'J which is quadratic *£ 

algorithm embodying this centering is calM RES1TC (suBP«Mly, « necessarily 

be‘en employed with good has never been 

JSS CJSSSaSo 4^ no. make any attempt . o treat the sutistics correctly, so 
that a it is not possible to assess the accuracy of the estimation adequately. 

A second approach has been put forth by Thompson [3,4], who 

sstessssssss: 

cessfully employed in support of the AMPTE spacecraft. 

Davenport [5] has proposed another 

i„g an app^ima^ produced by this 
the iterative solution of the least square p ^ i|s consistency. Hence, the ap- 

algorithm, unfortunately, mak “PP solution as the number of data becomes infinite, 

proximate solution cannot approach the true solution as tne n 

SHSHS'xrsc 

the magnetometers of the Hubble Space Telescope. 

by 

cohering; does no, suffer from th e ron vcrgence ^roWems o^a fb^-^myi^oridin^such^in 

Ref. 3 and 4; is much more direct than the algonth - ^ J U e completely and 

stage. It does this in several important ways by 1) treating t sta,.sn« m P ^ 

correctly (2) correcting for the centering operation, and ( 3 ) ^timating scale 
biases. The authors do not call this algorithm either RESIDA or RESIDS. 


The Model 

All treatments begin with the model ^ ^ 

where B fc is the 

tomt'^rtlbfcred coordinate system; -4 t is the < ^^^ m ?®JJ^ t ^^[ r ^jJ,ent S iwhe W 'The 

Earth-fixed coordinates; b is the “^“^erroriandgMmagnetic field model uncertain- 

.0 be whim and Gaussian. This is probably a poor approximation, 
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r ® rron5 in th u e geomagnetic field model are certainly highly correlated, and in fact 

Sere 

S ea“ C a,,i,Ude ' We ,ranSP05e ,CrmS «I“ a,i °" <» a "0 


\H k \ 2 = \A k H k \* = lB k -b-e k \> 

= |BJ 2 - 2B* • b + |b| 2 - 2(B t - b) -e* + |e fc | 2 . 

If we now define effective measurements and measurement noise according to 

z k = |B*| 2 - |H fc | 2 , 
v k — 2(B fc — b) • — |c fc | 2 , 

then we can write 

z k = 2B* b- |b| 2 + v k , k = l,...,N. v , 7 

P ° i , nt f °L the derivation of a » of the algorithms. (Note that in equations (3b) 
, d (4) ’ ,s va,ue about which ‘he measurement is linearized and therefore must be 
interpreted as the sample value of the measured magnetic field and not a random van”we.) 

Even with the assumption that the original measurement noise is white and Gaussian the 
fan iThaT 3511 ^" 16111 ^ wWte * nd Gaussian ‘ ^ming that e fc is white and Gauss 

€*~^(0, EJ, (5) 

and 


(2ab) 

(2c) 

(3a) 

(3b) 

( 4 ) 


£{**«!} = 0 for kjti, 


then 


( 6 ) 


Pk = E{ v k } = -tr (E fc ) , (7 a ) 

a\ S E{ V 2 } - £ = 4 (B fc - b) T E fc (B* - b) + 2 £(E fc )?. , (7b ) 

1 = 1 

so that v t must contain both Gaussian and x ! components, as is evident from equation (3b), 
Here tr( ■ ) denotes the trace operation. In addition, 

E i v k v e} = Vk»t t (8) 

so that the v* are uncorrelated but not white. If we assume, hoover, that the noise €fc is small 

a“pr P o« r maL°.y ge ° magnetiC *" * 3 ^ is Gaussian and - can Tte 


v k ~N{n k , <r 2 k ), 

keeping only the first term in equation (7b). 


( 9 ) 
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Maximum Likelihood Estimate of the Bias and Scoring 


Given the statistical model above, the negative-log-likelihood function [6] for the magne- 
tometer bias is given by 

* b + |b| 2 - fi k ) 2 + log<r£ + log 2trj , 


•w = ^E[^- 2: 

k= 1 k 


( 10 ) 


which is quartic in b. The maximum-likelihood estimate maximizes the likelihood of : th - 

mate which is the probability density of the measurements (evaluated at their sampled val ) 

given as a function of the magnetometer bias. Hence, it minimizes the negative logarithm of 
the likelihood (equation (10)), which thus provides a cost function. 

Since the domain of b has no boundaries, the maximum-likelihood estimate for b, which 
we denote by b*, which minimizes the negative-log-likelihood function, must satisfy 

= 0 . ( 11 ) 


dJ 

db 


Note that only the first of the three terms under the summation depends on the magnetometer 
“less L wishes .0 estimate parameters of .he measurement 

retain the remaining two terms. This quartic dependence can be avoided if complete three axis 
attitude information is available, since the bias term then enters linearly into the measurement 
model (q.v. equation (1)) as in the work of Lemer and Shuster [ 7 

The most direct solution is obtained by scoring, which in this case is the Newton-Raphson 
approximation. We consider the sequence 

»r=o, (12) 

This series is obtained by expanding J(b) to quadratic order in (b - b? R ), setting the gradient 
S^n^teS series to So, absolving for b i+ . If for some value of . we are suffictently 
close to the maximum -likelihood estimate, then it will be true that 

lim b^ R — » b* . ( 13 ) 

i— ► OO 

We have made the convention here that the partial derivative “5 “ 

to a column vector is again a column vector. The gradient vector dj/db is the 3 x 1 matrix 


d ± = _ y -L ( 2fc - 2B* • b + |b| 2 - /x fc )2(B fc - b) , (14) 

and the Hessian matrix d 2 J/dbd b T is given by the 3 x 3 matrix 


d2j - V — [4 (B fc - b)(B fc - b) T + 2 (z fc - 2 • b + |b| 2 -njhxt] ■ ( 15 ) 

dbdbT A, ff 2 l 

»■ Throughout thi work we shall use k as the time index and . as the iteration index. 
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Generally, the second term in the brackets will be much smaller than the first and can be 
discarded. 


A second approach to scoring is the Gauss-Newton approximation [8]. In this case, we 
replace the Hessian matrix by its expectation, the Fisher information matrix F. Since 

■£{(** T 2 B fc • b + |b| 2 — fi k ) } = 0 , (16) 

this amounts to discarding the second term. According to the law of large numbers, as the 
number of independent identically distributed (i.i.d.) samples of a random variable becomes 
infinite (the asymptotic limit), the average of these samples approaches the expectation value 
of the random variable. Our measurements are not identically distributed because of the de- 
pendence on B fc . However, if the distribution of the values of is regularly repeated, 
then we may regard the measurements as being i.i.d. for each value of A*H fc . Except for 
the replacement of the Hessian matrix in equations (12) by the Fisher information matrix, the 
iteration proceeds as before. 


For both the Newton— Raphson and the Gauss— Newton method, the estimate error covari- 
ance matrix is given in the limit of infinitely large data samples by 



E^( B *-b)(B,-b) T 

Jt=i k 


(17) 


If the measurement noise is Gaussian, then the asymptotic limit is true, in fact, for finite data 
samples. In most cases, the Fisher information matrix is simpler to evaluate than the Hessian 
matrix of the negative-log-likelihood function, and often can be evaluated independently of 
the data. 


The earliest estimates of the magnetometer bias were accomplished by the method culmi- 
nating in equations (12) though usually the weights were not chosen according to a statistical 
criterion. 


The Centered Estimate 


In order to avoid the minimization of a quartic cost function, let us define the following 
weighted averages 



This is similar to the centering approximation of Gambhir [1,2], who, however, did not deter- 
mine the weights from any statistical quantities. It follows that 


z = 2Bb — |b| 2 + v. (20) 
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If we define now 


z k = z k — z , B fc = B*-B, v k =zv k -v, ji k = n k fi, (21) 

then subtracting equation (20) from equation (4) leads to 

z fc = 2 B fc • b + , k = l,...,N. (22) 

This operation is called centering. 

The centered measurement is no longer quadratic in the magnetometer bias vector. How- 
ever, the centered measurement noise is no longer uncorrelated. Thus, one can no longer write 
the negative-log-likelihood function in the form of equation (10). Nonetheless, in practice one 
has ignored this and determined the bias from a cost ftinction of the form 2 


japprcK (b) = \ £ X (i k ~ 2B k • b - fi k ) 2 , (23) 

1 Jt=l * 

and achieved reasonable results in spite of the lack of consistency, arguing that one was only 
discarding a single measurement out of many. We shall see below that one can discard much 
more than 1 /N of the accuracy by this operation, but we shall see also that equation (23) is 
closer to being correct than one might have imagined. Note that the sum is from 1 to N - 1, 
since the centered measurements are not independent. 

£> = °. (24) 

k=l 


Minimizing t/ approx (b) over b leads to 

= p;r~ E A ( 5 * - , (25) 

fc=l k 

with the estimate error covariance matrix given approximately by 



Note that fi k will vanish if the original measurement noise e k , k = l, , N, is identically 
distributed. The centered estimator converges in a single iteration because the cost function 
is exactly quadratic. 

*111 actual feet, these calculations have almost always assumed a constant weighting and neglected the contri- 
bution of nt- 
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Fixed-Point Method 


1b avoid the loss of data from centering, Thompson, Neal and Shuster [3,4] proposed a 
fixed-point algorithm. Define the quantities 


G = [ 4B * B I + 2 (*fc-/**H3x3] , (27a) 

k = 1 * 

N 1 

a = ’( z fc-^)2B fc , (27b) 

*=i k 
N 1 

f ( b > s E;r [4(B t .b)b + 2|b|*(B,-b)] . (27c) 

*=1 * 

Then the gradient of the negative-log-likelihood function becomes 

dJ( b) 

-^- = Cb-a-f(b) = 0. (28) 

which can be solved implicitly to yield 

b* = G~ l [a -f- f(b*) ] . (29) 


This equation must be solved iteratively, 

l>o FP = 0, b&-Gr->[« + f(bP)], (30) 

and we expect that once bf 1 * is sufficiently close to the solution that 

HmbP = b\ ( 31 ) 


Davenport’s Approximation 


Davenport and his collaborators [ 5 ] have offered an approximate form for the bias vector 
estimator. He begins by writing an approximate cost function as 

1 N 1 

^(b) = 2 £ ^2 ( z k ~ 2 B fc • b + A 2 - fi k ) 2 , (32) 

z fc=i a k 

where A is a constant. This cost function would agree with that of equation (10) when A = |b|. 
Davenport, however, allows A to be a free parameter. 
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The cost function of equation (32) is only quadratic in b. Differentiating this cost function 
with respect to b and setting the gradient equal to zero leads to a solution of the form 

bJ> = f D (A), ( 33 ) 

that is, the estimate of the bias is a function of the parameter A. The “consistent value of this 
parameter is obtained by solving 

|fp(A)| 2 = A 2 . (34) 

Because Davenport’s algorithm effectively changes the dependence of the non-random part 
of the measurement on the bias even in the absence of noise, it cannot be consistent. Thus, 
as more data is accumulated the accuracy will not improve. However, it can be used as the 
starting point for a Newton-Raphson or Gauss-Newton iteration of the quartic cost function. 


A Statistically Correct Centered Algorithm 

The original data, z k , k = 1, ... , N, may be replaced by the centered data, z k , k = 
1 N - 1, and the center value z. The measurement equations are given by equations (20) 
and (22). The centered data have the advantage of depending only linearly on the magne- 
tometer bias. However, they have the disadvantage that the centered measurement noise is 
correlated. Therefore, the cost function for the centered data alone cannot be written as the 
sum of N - 1 squares. To write a statistically correct cost function foUhe centered data (making 
the approximation that the measurement noise v k is Gaussian) we define 


Z = [ , Z 2 , . • 

i T 

• > Z N- 1 J » 

^ ^ ^ T 

B = [B 1? B 2 , ... , , 

(35ab) 

■M S [fill P"2i • 

1 T 

• ■ 7 Pn- 1 J > 

V = 1 ] » 

(35cd) 

and write formally 

2 = 2Bb + V, 

(36) 

with 

V ~ jV f-M, • 

(37) 


Here % is the covariance matrix of V. (Note that B is an (N - 1) x 3 matrix) 

The negative-log-likelihood function for this stacked centered measurement is simply 

J (b ) = I \(z-2Bb-M) T it- 1 (z - 2Bb - A4 ) + log det 72 + (N - 1) log 2x j . 

2 L V 1 (38) 

Equation (23) made the assumption that U was diagonal. We do not make this approximation 
here. Minimizing this negative-log-likelihood function leads directly to 

b* = (4 6 r ^- 1 5) _1 2B v U~ l (Z-M), (39) 
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with estimate error covariance matrix 


p b b = (iiFk-'b) 1 


(40) 


quantities of data, the naive evaluation of equations (39) and (36) can be a for 
S“£ ‘ aSt T1,erCf ° rc ’ we “* ,|K — inverting the matrix in U«on ’ro £u* 

K kt = E{ (v k - ji k )(v t -£,)} = a\ 6 k( - a 2 , ( 41 ) 

which shows the correlation explicitly. However, this matrix has the simple inverse 

(*")„= ^ (42) 

where a% is the variance of v N . Substituting this expression into equation (38) leads to 
1*1 

= 2 E 7^*" 2 ®* ' b ~ £*) 2 + terms independent of b . 


k = 1 v k 


(43) 


The statistically correct cost function for the centered data looks exactly like the naive exnres 

‘ ha ‘ ' he SUmmMi0n * "° W fr ° m 1 “ N - The minimi “"'on is 


~ ~ 1 

^ ~ ^ 53 (*k ~ fik) 2B k , 

k = 1 a k 

and the estimate error covariance of the centered estimate is given by 



(44) 


(45) 


anJk tCred eS !‘ rT ! ate 1S now to ** much m ore attractive than before. It is simple 
and by a very trivial alteration (replacing the sum from 1 to N - 1 by a sum from 1 to N) it 

can be made to treat the statistics of the measurement noise correctly. It is very different 
n character from that the centered estimate of Gambhir [1,2]. It is thus to be preS 
to^ompson s algonthm [3, 4] , whose convergence can be problematic, and to Davenport’s 

IS T [5 ] ’ W r h “ n0t consistenL ^ greatest drawback to the centered algorithm 
lies m the exclusion of certain data, the effect of which we now investigate 8 


The Complete Solution 


nr *°* N la ' ge, 3 h ? nai ? centered a, gorithm presented earlier is hardly worse than the rig- 
usly centered algorithm derived above. From the standpoint of computation burden, the 
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mote rigorous treatment of the statistic* has merely added one more term (out of W) to the 
summation. However, equation (45), because it has been derived rigorously, will afford us the 
possibility of computing the correction from the discarded measurement z. 

Instead of the measurement set z k , k = 1, . .. , iV - 1, 7, we may now consider ^emea 
surements to be effectively b* and z. Therefore, to determine the exact maximum likelihood 
estimate b*, we must develop the statistics of these two effective measurements more com- 
pletely. 

Let us substitute equation (22) into equation (44). This leads to 


N 


b* = E 4(2B fc -b + t> fc -Mfc)2B fc . 


(46) 


— at 
k=i k 


which we may rewrite as 


jv 


b* = b + P bb 53 i 2Bfc (®k - H) 

k = i ** 


= b + v fc . 

The last term is just the (zero-mean) estimate error. Obviously, 


v 6 ~ A r (o, P bb ) 


It follows that we can write 


(47a) 

(47b) 

(48) 

(49) 


j(b) = I (b - b*) T P bb (b - b*) + terms independent of b , 

which can be verified by expanding equation (43) and completing the square in b. The estimate 
b* is thus a sufficient statistic for b [6] . Equation (49) is very useful, because it allows us to 
investigate the effect of corrections to the centered formula using only our knowledge of b 
and P. We do not have to refer again to the N centered measurements z k , > ■ • • » 

We must now combine b* and J to obtain a complete representation of our data for the 
computation of b. Recall equation (20), 

( 20 ) 


z = 2B • b - |b| 2 + «, 


with 


Note that z which unfortunately, is a nonlinear function of b, is nonetheless an extremely 
« m^mfment, more accurate than the other 

1 UfN because a is smaller typically than the other variances by this factor. Thus, simply 
ca'nuil asignificant loss of accuracy if 5 - b is not significantly smaller 

than typical values of B fc . 
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(51a) 


What is the correlation between and v? Calculating this explicitly, gives 

N 1 

£{ v,(» - » ) = p Y, ^ B*«{ («* - W(® - » } 

*=x * 

-^1- 

= P E ^V = 0, (51bc) 

Jt=i * 

and we have used equation (21). Thus, v t and v are uncorrelated. It follows, that the negative- 
log-likelihood functions add and 

7(b) = 7(b) + 7(b), (52) 

with J (b) given by equation (49) and 

*^(b) = 2 = 2 ( 2 - 2B b + |b| 2 -ft) 2 +logff 2 + log2xJ (53) 

The weight associated with the center term 7(b) is equal to the sum of all the weights of 7(b). 
Thus, when B is not small, the loss of accuracy from discarding the center time can be sub- 
stantial. We can determine the relative importance of these terms to the estimate accuracy by 
computing the Fisher information matrix F bb to obtain 

Fh = Ph + £ (5 - b)(B - b) T = P? . (54) 

The estimate error covariance matrix will be the inverse of this quantity. If the distribution of 
the magnetometer measurements is “isotropic,” that is, if B - b vanishes, then 7(b) will be 
insensitive to b. It is in this case that the centering approximation obviously leads to the best 
results. If, however, one attempts to determine the magnetometer bias from a short data span, 
say, from an inertially stabilized or Earth-pointing spacecraft, then B - b will be equal to the 
similar expression for a typical value of the magnetic field, and the formerly discarded center 
term which will provide half of the accuracy, especially for the component along B — b. 

Thus, our new algorithm is as follows: 

• We compute the centered estimate b* of the magnetometer bias and the covariance matrix 
P bb using the centered data and equations (44) and (45). 

• Using the centered estimate b* as an initial estimate, the correction due to the center term 
is computed using the Gauss-Newton method 

b.+^b.-F^Cb,)^), (55) 

where the Fisher information matrix F bb ( b) is given by equation (54), and the gradient 
vector is given by the sum of the gradients of equations (49) and (53) 

g(b) = (b - b*) - ± (z - 2 B • b + |b| 2 — Ji) 2 (B - b) . (56) 
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• The iteration is continued until 

Vi = (b< " b i _ 1 ) T n6(b i -i)(b j - b^j) (57) 

is less than some predetermined small quantity. 


Numerical Examples 

The algorithms treated in this work have been examined for an inertially stabilized space- 
craft The spacecraft orbit has been chosen to be circular with an altitude of 560 km and an 
inclination of 38 deg. This is, in fact, the orbit of the SAC-B spacecraft (Satelite de Aplica- 
ciones Cientificas), the first spacecraft to be developed by Argentina, which will be inertially 
stabilized to observe the Sun. The geomagnetic field in our studies has been simulated using 
the International Geomagnetic Reference Field (IGRF (1985)) [9], which has been extrap- 
olated to 1994. More recent field models are available, but IGRF (1985) is adequate for our 

simulation needs. 

for purposes of simulation we have assumed an effective white Gaussian magnetometer 
measurement error with a standard deviation per axis of 2.0 mG, corresponding to an angular 
error of approximately 0.5 deg at the equator. We have assumed that no axis of the magne- 
tometer is predominantly parallel to the spacecraft spin axis or the geomagnetic field. The 
data were sampled once every ten seconds. 

We examine first Davenport’s approximation. To highlight the inconsistency of this method, 
we examine its behavior and that of the centered estimate for noise-free data. "Hie results for 
half an orbit of data for the spinning spacecraft are shown in Table 1. The equivalent results 
for noisy data are presented in Table 2. 

For small values of the magnetometer bias, Davenport’s approximation yields acceptable 
results. For values of the magnetometer bias comparable to or greater than the magnitude 
of the ambient magnetic field, the errors in Davenport’s approximation become unacceptably 
large. These statements hold both for the noise-free and the noisy data. 

We can gain a greater appreciation of the behavior of these two algorithms if we examine 
the normalized errors, f/ and tjd defined by 

Tbble 1. Comparison of Davenport’s Approximation and Centered Estimate for Noise-Free 


Data 


Model Bias (mG) 

Centered Estimate 

Davenport’s Approximation 

[ 10., 20., 30.] 
[ 30., 60., 90.] 
[ 60., 129., 180.] 
[100., 200., 300.] 
[200., 400., 600.] 

[ 10., 20., 30.] 
[ 30., 60., 90.] 
[ 60., 129., 180.] 
[100., 200., 300.] 
[200., 400., 600.] 

[ 10., 20., 30.] 
[ 30., 50., 90.] 
[ 60., 129., 180.] 
[101., 208., 261.] 
[180., 539., 161.] 
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Tbble 2 . Comparison of Davenport’s Approximation and Centered Estimate for Noisy Data 


Model Bias (mG) 

Centered Estimate 

Davenport’s Approximation 

[ TO., 20., 30.] 
[ 30., 60., 90.] 
[ 60., 129., 180.] 
[100., 200., 300.] 
[200., 400., 600.] 

[ 9.88 ±0.35, 20.47 ±0.82, 27.77 ±2, 44] 
[ 29.69 ± 0.30, 60.52 ± 0.69, 89.95 ± 2.03] 
[ 59.47 ± 0.25, 130.98 ± 0.65, 174.54 ± 1.58] 
[100.33 ± 0.22, 201.13 ± 0.56, 296.78 ± 1.60] 
[199.82 ± 0.37, 400.76 ± 6.34, 598.28 ± 3.44] 

[ 10.27, 19.88, 29.70] 
[ 29.66, 60.60, 89.66] 
[ 54.88, 129.28, 179.75] 
[101.94, 214.30, 230.13] 
[178,92, 538.68, 160.33] 


1 s [(•>"“ - b*) 1 ^ 1 <b— - b) - 3] , (58) 

which should have mean zero and standard deviation unity. A comparison of these quantities 
IS given in Table 3. The inconsistency of Davenport’s algorithm is evident. 

To see the advantages of this algorithm over naive quartic scoring consider the estimation 
of a magnetometer bias whose true value is (10., 20., 30.) mG. The results of successive it- 
erations for naive quartic scoring and the new algorithm are shown in Table 4. For the new 
algorithm, the first algorithm is the statistically correct centered algorithm and successive it- 
eration are obtained by applying the Gauss-Newton method to the complete cost function as 
given by equation (53) The l<r error brackets, computed from the Fisher information matrix, 
are (±.13, ±.19, ±.12) mG. The results of the two methods are nearly identical in this case 
and the convergence is equally rapid. Small differences in the results are due to the slightly 
different treatment of the roundoff errors. 

Consider now the case where the magnetometer bias vector is large compared with the 
ambient field, say (100., 200., 300.) mG. In this case we obtain the value presented in Table 5. 
Ine la error brackets here are found to be (±.12, ±.10, ±.12) mG. 

m this case naive quartic scoring does not even converge to the correct answer, nor does the 
method of Thompson, Neal and Shuster [3,4], which does not converge at all. Naive quar- 
tic scoring converges, in fact, to a local minimum. The new algorithm, on the other hand, works 


Thble 3. Comparison of Normalized Errors for Davenport’s Approximation and the Centered 
Estimate for Noisy Data 


Model Bias (mG) 

n 

Vd 

[ 10., 20., 30.] 

0,82 

- 1.18 

[ 30., 60., 90.] 

0.68 

3.87 

[ 60., 129., 180.] 

0.814 

0.81 

[100., 200., 300.] 

2.66 

2.86 xlO 3 

[200., 400., 600.] 

-0.83 

0.65 xlO 3 
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•ftbl. 4. Comparison ofN.lv* Quartlc Scoring and the New Algorithm. Tbs true value of the 
magnetometer bias vector Is ( 10 ., 20 ., 30 .) mG. 


Iteration 


Naive Quartic Scoring 


1 

2 

3 

4 


[ 10 . 08 , 19 . 27 , 33 . 04 ] 
[ 9 . 84 , 20 . 18 , 29 . 91 ] 
[ 9.84 , 20 . 19 , 29 . 89 ] 
[ 9.84 , 20 . 19 , 29 . 89 ] 


New Algorithm 


[ 9 . 82 , 20 . 08 , 29 . 05 ] 
[ 9 . 90 , 19 . 83 , 29 . 94 ] 
[ 9 . 90 , 19 . 83 , 29 . 93 ] 
[ 9 . 90 , 19 . 83 , 29 . 93 ] 


vety well. Note that a single iteration of the center correction is sufficient. Hie errors for the 
new algorithm are clearly consistent with the computed confidence intervals. 


Discussion 

The new algorithm for attitude-independent magnetometer bias determination produc^ 
excdlentTesults in all situations. Since it begins with a very good initial estimate for the 
bias, it is more likely to converge to the correct minimum than does oa.vesco-'ngt 1J cn the 
fixed-point method of Thompson et al [ 3 , 4 ], which begin at b = 0. Unlike the centered 
i nf RFSIDG fame r 2 1 it does not discard data and does the centering m a 

Sw^cct waf^. fiom S approbation tha, the measurement errors on the attitude- 
independent derived measurement are Gaussian and uncorrelated, which is almost «rtaiidy 
not the case It is amusing to speculate that the statistically correctly centered cost function 
of equation (43) would probably be rejected as statistically incorrect by heuristic algorithm 
devdoperT unschooled iS Statistics, because it appears to use redundant data. Its irntial cen 
tered Jstimate for the magnetometer bias is clearly a better approximation than ignoring the 
aua1r“ha^or of Urn in the work of Davenport et al [5]. The new algorithm is cer- 
tain! more sophisticated statistically than its predecessors, and more efficient compu ationa ly^ 
Perhaps, most importantly, the new algorithm makes manifest the physical quantities which 

Thble 5. Comparison of Naive Quartic Scoring and the New Algorithm. The true value of the 
magnetometer bias vector is ( 100 ., 200 ., 300 .) mG. 


Iteration 

1 
2 

3 

4 

5 

6 
7 


Naive Quartic Scoring 


[ 107 . 62 , 259 . 77 , 
[ 51 . 51 , 398 . 62 , 
[ 70 . 35 , 358 . 17 , 
[ 72 . 13 , 340 . 88 , 
[ 71 . 78 , 338 . 71 , 
[ 71 . 70 , 338 . 64 , 
[ 71 . 70 , 338 . 64 , 


2 . 85 ] 
- 368 . 88 ] 
- 196 . 33 ] 
- 145 . 65 ] 
- 140 . 60 ] 
- 140 . 62 ] 
- 140 . 62 ] 


New Algorithm 

[ 99 . 82 , 200 . 63 , 298 . 02 ] 
[ 99 . 97 , 200 . 11 , 299 . 81 ] 
[ 99 . 97 , 200 . 11 , 299 . 81 ] 
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determine the behavior of the bias estimator. It is also understandable. 
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